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Due to the large tt production cross-section at the LHC energies, the ATLAS experiment is expected to have 
enough statistics to measure tt cross-section even at initial luminosities. Recent studies performed in ATLAS 
on the development of tt cross-section measurements in the lepton+jets channel at y/s=lO TeV using the 
likelihood method will be discussed. The expected statistical and systematic uncertainties for the cross-section 
measurement using the likelihood method are evaluated for an integrated luminosity of 50 pb~^ of Monte Carlo 
(MC) simulated data. Measurements with data that will be collected in the first year of the LHC operation are 
emphasized. 



1. Introduction 

The top quark was discovered in 1995 at Fermilab in 
the pair production mode (tt events) through strong 
interactions [1]. At the LHC, production takes place 
at small momentum fractions of the proton, x w 10~^ 
where the gluon parton distribution function (PDF) is 
large. Therefore, the gluon scattering process (« 90% 
at 10 TeV) dominates over quark-antiquark scattering 
in tt production at the LHC. Because of much larger 
parton luminosities (dominated by the gluon PDF), 
the tt production cross-section, o-jf, at the LHC is 
significantly larger than at the Tevatron. Assuming 
mt = 172.5 GeV and using CTEQ6.6 and CTEQ6.5 
PDFs [2], the ti cross-section at 10 TeV is 400 pb±ll% 
at next-to-leading order (NLO) and 400 pb±6% at 
next-to-next-to- leading order (NNLO). Any deviation 
in the measured cross-section from the theoretical 
value can be considered an indication of the presence 
of physics beyond the Standard Model. Measurements 
of the tt cross-section in different channels (e-|-jets, 
fi+jets, T-|-jets, di-leptons) provide a more stringent 
test of Standard Model predictions, since new physics 
may modify relative fractions of these channels with 
respect to the total cross-section. Top quark events 
will be a major background for Higgs and supersym- 
metry (SUSY) searches. Therefore, accurate knowl- 
edge of the tt cross-section is crucial during the early 
operation of the LHC to discover new physics. 

In the Standard Model, the top quark decays almost 
exclusively through the t — + Wb decay. From the ex- 
perimental point of view, we consider the case when 
one of the W^'s decays leptonically, and the other de- 
cays hadronically (tt WWbb l^'jij2bb with I = 
e, /i and ji, j2 are light jets). This is known as a semi- 
leptonic tt decay. This allows us to trigger events 
efficiently using the single lepton trigger. The final 
state of the semi-leptonic channel consists of a high- 
Pt charged lepton, missing energy and four jets. This 
decay channel is usually called lepton-|-jets, and we 
will use this terminology in this paper. 

Several methods for the ti inclusive production 



cross-section measurement have been developed and 
carefully studied by the ATLAS collaboration [3] . All 
of these methods have the goal of an early tt cross- 
section measurement. Measuring the tt cross-section 
using the likelihood method can be complementary 
to other analyses and therefore provides a valuable 
cross check of the early results. This method uses 
shapes of different kinematic variable distributions in 
ti and background events. We choose kinematic vari- 
ables which have good separation power between ti 
and backgrounds. We optimize those discriminating 
variables to reduce statistical and systematic uncer- 
tainties. We have not used b-tagging [3] in this anal- 
ysis. 

2. The Likelihood Discriminant 

The discriminating function can be generally writ- 
ten as follows: 

^ ^ S{X1,X2, ...,XAf) 

S{X1,X2, ■■■,Xn) + B{xi,X2, ■■■,xn) 

where xn are discriminating variables and S and B 
are N-dimensional probability density functions for 
signal and background, respectively. We can assume 
to good approximation that the N discriminant vari- 
ables are uncorrelatcd. Therefore, each multidimen- 
sional probability density function can be expressed 
as a product of one dimensional probability density 
functions: 

r _ ^iLlSiiXj) 

exp(Et,ln|fe)) (2) 

exp(Ef=iln|^) + l 

Denoting the fit to the logarithm of S and B as 
(In ■^){***'^''(a;i), we can write the likelihood discrimi- 
nant as 

^ exp(Eti(ln|)r"'^(xO) 
exp(Eti(ln|)r"^(xO) + l 



2 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



3. Cross-section extraction procedure 



Jets 



The U production cross-section is given by 

Ntt 



Br ■ L - Csei 



(4) 



where 7Vj( is the number of selected tt events, Br is 
the branching ratio for the final state considered, L 
is the integrated luminosity and Cgei is the selection 
efficiency. 



4. Reconstruction of physics objects 

We have used semi-leptonic decays of ti events 
(tt — > lz/jij2bb with I = e,ii) for this analysis. In this 
section, we define basic reconstructed objects (elec- 
trons, muons, jets and missing transverse energy) used 
in this analysis. These definitions are based upon 
the standard reconstruction algorithms recommended 
in ATLAS. The pseudorapidity is defined as rj — -In 
tan(|) where 9 is the polar angle from the beam axis. 
The azimuthal angle cj) is measured around the beam 
axis. We define separation of two objects m rj — <f) 
space as Aii = -/(A^^qiA^. 



Electrons 

Electrons are reconstructed in the inner tracker and 

in the calorimeter of the ATLAS detector. Electrons 
are required to have pseudo-rapidity in the range \ri\ < 
2.47 and px > 20 GeV. Furthermore, if an electron is 
foimd in the calorimeter crack region 1.37 < \ri\ < 
1.52, the electron is discarded. The electron has to 
be isolated based on the calorimeter energy: the ad- 
ditional transverse energy {Et) in a cone with radius 
AR = 0.2 around the electron axis is required to be 
less than 6 GeV. 



l\/luons 

Muons are reconstructed by combining information 
from the inner detector and the muon spectrometer. 
The muons used in this analysis must have minimum 
transverse momentum px > 20 GeV and are required 
to lie in the pseudo-rapidity range \r]\ < 2.5. Muons 
are isolated based on the calorimeter energy: the ad- 
ditional transverse energy Et in a cone with radius 
AR = 0.2 around the muon is required to be less 
than 6 GeV. In order to remove muons coming from 
decays of hadrons inside jets (i.e., i?-hadrons origi- 
nating from the ?>-quark in top quark decays), muons 
which are within a cone of size AR < 0.3 from a jet 
are removed. 



Jets are reconstructed using the standard ATLAS 
cone algorithm [3] in r] — ^ space, with a cone radius 
of 0.4. The jets used in this analysis are required to 
have minimum transverse momentum pt > 20 GeV 
and are required to lie in the pseudo-rapidity range 
\ri\ < 2.5. Jets which overlap with electrons within a 
cone of size AR < 0.2 are discarded. 

I\1issing Transverse Energy 0t) 

For the calculation of missing transverse energy 
we determine the sum of the following components: 
the contribution of cells in identified electron or pho- 
ton clusters, the contribution of cells inside jets, the 
contribution of cells in topological clusters outside 
identified objects, the contribution from muons and 
the cryostat correction [3]. The sum of total trans- 
verse energy in semi-leptonic top quark events is about 
500 GeV, which gives a typical resolution of the 
order of 10 GeV [3]. 



5. Selection of lepton + jets events 

We identify semi-leptonic tt events by requiring that 
either the single isolated electron or single isolated 
muon trigger, both with 15 GeV thresholds, have 
fired. The trigger efficiencies can be measured either 
from MC simulation or from data (Z — > ee or Z ^ /z/x) 

Apart from the trigger requirement, the following 
additional event selection cuts have been applied for 
event selection: 

• Exactly one isolated lepton (electron or muon) 
with Pt > 20 GeV. 

• ^T >20 GeV. 

• at least 3 jets with pt >40 GeV. 

• at least 4 jets with pt >20 GeV. 

The fraction of all tt events passing the individual 
selection requirements and overall commulative effi- 
ciency are shown in Table I. 

Table 1 Cut flow efficiency for the tt decay. 





/i+jcts (%) 


c+jots (%) 


Trigger 
Isolated lepton 

> 20GcV 

3 jets with pr > 40 GeV 

4 jets with pr > 20 GeV 


31.0 ±0.1 
67.0 ±0.1 
91.2 ± 0.1 
47.7 ± 0.2 
80.0 ± 0.2 


25.0 ±0.1 
70.3 ± 0.2 

90.6 ±0.1 
47.5 ± 0.2 

79.7 ± 0.2 


Cumulative efficiency 


7.23±0.04 


6.01±0.04 
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6. Backgrounds 

The lepton+jets data have a mixture of both ti 
events and background events (T-V+jcts, multi-jet 
QCD, single top and di-boson). The expected num- 
bers of tt and backgrounds in 50 pb~^ are shown in 
Table II. The major background is W+jets. In order 
to construct the likelihood template, we used discrim- 
inating variables which have good separation power 
between the tt and PF+jets events. The di-boson 
background behaves like the W+jets background, and 
therefore it is included as part of the VK+jets tem- 
plate. The single top likelihood template behaves 
like the tt template. Therefore, the number of tt-like 
events obtained with the likelihood templates fit in- 
cludes the single top production. We will subtract the 
expected number of single top events from the signal- 
like events using the theoretical cross-section of single 
top production. QCD backgrounds are expected to be 
present in the data sample due to imperfect jet and 
lepton identification and due to the presence of heavy 
flavor jets, where & or c quarks decay semi-leptonically. 
The number of expected QCD events will be deter- 
mined solely from data using the matrix method [5]. 
The contribution of QCD events will be subtracted 
from the data in each bin of the likelihood distribu- 
tion. 

Table II Expected numbers of ti and major back- 
ground events in 50 pb~^ of Monte Carlo (MC) 
simulated data. 



Data sample 


/i+jots 


f:'+jctS 


tt 


78.") 


()53 


VK(/i!/)+jots 


416 





W{ep)+iets 





297 


single top 


56 


52 


Di-boson 


3 


2 


QCD: Di-jet 


8 


13 



7. Likelihood Discriminant variables 

For the likelihood function, we choose variables that 
satisfy the following requirements: 

1. Their distributions in data arc expected to be 
accurately reproduced by MC simulation; 

2. They have either minimal or no dependence on 
the jet energy scale corrections which could be 
a source of large systematic errors in early data; 

3. They have either minimal or no dependence on 
missing energy; 

4. They are not highly correlated with each other. 



The following topological variables have been used to 
construct the likelihood function : 

• Centrality = ; Ht is the scalar sum of the 
Pt of the four leading jets and H is the scalar 
sum of the energy of the four leading jets, 

• Aplanarity A = | x A3 ; A3 is the smallest 
eigen value of the normalized momentum tensor 

Mij = , where p° is the momentum vec- 

^ Eol(p°)P 
tor of reconstructed object o (e, /i and four lead- 
ing jets) and i and j are Cartesian coordinates. 
After proper diagonalization, Mij can have three 
eigenvalues, Ai > A2 > A3, with A1-I-A2-I-A3 = 1. 
Aplanarity is the measure of the flatness of 
the events. Spherical events should have large 
A whereas planar events should have small A. 
W+jets events are more planar (small A) than 
ti events. 

• Sphericity S' = | x (A2 -I- A3); where A2 and A3 
are the smallest eigenvalues of the normalized 
momentum tensor (A3 < A2) so that < 5 < 1. 
Sphericity is a measure of the summed with 
respect to the event axis. For isotropic events, 
5 Ri 1 and for two jet events 5 ~ 0. 

• r]e where i = e,n; 

• A77(j2, J3) = vij2) - viis); 

• Fi9(;,9j,) = 2 I tan-i(e-''W)-tan-i(e''^^'^') |; 

• FiOi^Oj^) = 2 I tan-i(e""*^')-tan-i(e''^^''') |; 

• F{9e,ej,) = 2 I tan-i(e-''(^))-tan-i(e"(j'^^)) |; 

The likelihood templates obtained using these topo- 
logical variables are shown in Figure 1 and Figure 2. 

8. Expected statistical uncertainty 

In order to estimate the statistical uncertainty, we 
performed several ensemble tests where we mixed dif- 
ferent proportions of tt and VF-|-jets events to pro- 
duce "pseudo data" . We would like to see whether we 
can extract the exact fraction of tt events from the 
"pseudo data" using the default likelihood template. 
We have divided the tt samples into two halves. The 
first half of the tt samples was used for likelihood tem- 
plate construction. The second half of the ti sample 
was used for the ensemble tests. Due to the limited 
statistics of the VF-l-jets samples, we used all available 
VF-hjets samples for both VF-hjets likelihood template 
construction and the ensemble test. We performed 
1000 trials to obtain a distribution of the fitted signal 
fraction. Figure 3 and Figure 4 show the ensemble 
fit for the muon channel and electron channel respec- 
tively. 
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Figure 1: Likelihood templates for ^+ jet channel: 
template (solid red line), W+jets template (blue 
dashed line) 









0.22 




ATLAS Work in Progress 




0.2 






f 0.18 
1 








|0.16 
















|0.14 








0.12 






0.1 














0.08 




: 








:, , , 


. . . . 1 .... 1 . . . . 1 .... 1 .... 1 .... i .... 1 .... 1 ... . 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 2: Likelihood templates for e+jet channel : 
ti template (solid red line), VK+jets template (blue 
dashed line) 
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The initial tt fraction, average fitted fraction of tt 
events from the ensemble fit, and RMS of the dis- 
tribution are presented in Table III. There is good 
agreement between the input and fitted values of the 
tt fractions. The RMS of the distribution of the fitted 
fraction is considered to be the expected statistical 
error of the method. 



Table III Initial tt fraction, average tt fitted fraction 
and RMS for ensemble test 



Decay Channel 


Initial 
fraction (%) 


Fitted 
fraction (%) 


Statistical 
error (%) 


e+jet 


69 


68 


6.0 


M+jet 


65 


64 


5.3 



9. Systematic uncertainties 

The systematic uncertainties on the tt production 
cross-section (Acrtj) originate from uncertainties on 
the selection efficiency (Aesej) and uncertainties on 
Nti {ANtt, from the variations of the likelihood tem- 
plates in the fit). Considering the fact that the sys- 
tematic uncertainty from a given source can affect 
both the selection efficiency and the shape of the like- 
lihood template at the same time, the systematic un- 
certainty on the tt production cross-section is deter- 
mined by varying the source by one standard deviation 
up and down and propagating the variation into both 
fitted number of tt events and the signal efficiency [5] : 



final state radiation (FSR), and PDF are discussed 
below. 

9.1. Jet Energy Scale (JES) 

To obtain the uncertainty due to the jet energy 
scale, we varied all jet energies coherently by 10% 
(pessimistic senario) and 5% (optimistic senario). 
has been corrected after varying the jet energy scale. 
We have seen that the likelihood template does not 
change significantly when we vary the jet energy scale. 
We calculated the modified selection efficiency for tt 
events after applying the jet energy scale. The change 
in selection efficiency due to variation of the jet energy 
scale propagates as a systematics unccrtainity to the 
cross-section. 

To estimate the systematic uncertainty due to the 
likelihood shape, we derive new likelihood templates 
for tt and VF-|-jets after changing the jet energy. We 
make "pseudo data" using samples after varying the 
jet energies. Then we use the nominal likelihood tem- 
plate to find out the exact fraction of tt events from 
the "pseudo data" . The difference between fitted frac- 
tion of ti events and the original fractions of tt events 
is considered to be the error due to the likelihood 
shape for JES. The systematic error due to the like- 
lihood shapes for JES is negligible compared to the 
error due to selection efficiency. The contribution of 
JES systematics to the has been summarized in 
Table IV. 

9.2. MC generators 



ftt ± Aa-tt — (5) 

Br ■ L ■ [esei ± l^esei) 

Major sources of uncertainties due to jet energy 
scale, the Monte Carlo model used for the signal and 
background simulation, initial state radiation (ISR), 



We calculated the uncertainty due to the MC gen- 
erators in a similar manner as the JES systematic un- 
certainty. We used ACER MC to make a new "pseudo 
data" and applied "nominal" templates to this data. 
The difference between the fitted fraction obtained on 
ACER and MC@NLO events is considered to be the 
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Probability 



MC generator systematic uncertainty. We also calcu- 
lated the modified selection efficiency for ti events and 
used it as the MC generator systematic uncertainty 
on the selection efficiency. The total systematic un- 
certainty due to MC generators is shown in Table IV. 



9.3. ISR/FSR Modeling 



More initial state radiation (ISR) and final state 
radiation (FSR) increases the number of jets and af- 
fects the transverse momentum of particles, tt selec- 
tion cuts include these quantities. Therefore, ISR and 
FSR will have an effect on the selection efficiency. In 
order to estimate the effect of the ISR/FSR systemat- 
ics on the selection efficiency, we compared the results 
obtained with the tt AcerMC-generator [6] interfaced 
to PYTHIA [7] for showering and hadronisation. The 
Pythia parameters which control the ISR/FSR were 
then varied in order to generate the maximum dif- 
ference for the reconstructed hadronic top masses and 
two reference samples created. The relative differences 
between the numbers of selected events obtained with 
these reference samples, relative to the default central 
sample, was taken as the uncertainty on the cross- 
section. In order to estimate the effect of ISR/FSR 
systematics on the likelihood shape, the nominal like- 
lihood templates are obtained from half of the default 
central sample. Using the second half of the cen- 
tral sample and W^-l-jets backgrounds, we produced 
"pseudo data", then we used the nominal likelihood 
template to find tt fractions from the "pseudo data" . 
A similar procedure was applied to two reference sam- 
ples to find tt fractions from the "pseudo data" and 
the difference between the mean of the nominal sample 
to the mean of the reference samples was considered 
to be the uncertainty due to ISR/FSR. The system- 
atic error due to ISR/FSR has been summarized in 
Table IV. 




^0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Probability 



9.4. Parton Distribution Function (PDF) 
uncertainties 

The uncertainties due to the PDF are examined by 
a re-weighting scheme that uses MC truth informa- 
tion about the hard partons. The basic procedure is 
to evaluate the probability of an event with a partic- 
ular kinematic characteristic to be produced [4] . The 
variation of the selection efficiency on the signal sam- 
ple was taken as a measure of the cross-section mea- 
surement variation [8]. The errors coming from the 
PDF uncertainties are shown in Table IV. We did not 
evaluate the PDF uncertainty due to the change of 
the likelihood shape since this error is significantly 
smaller compared to others, like JES, ISR/FSR and 
Monte Carlo generators. 

9.5. Trigger efficiency 

The lepton trigger efficiency is measured from the 
simulated MC data using Z events. We expect the 
uncertainty to be of the order of 1% for electrons and 
muons [4]. 

9.6. Single top 

The systematic uncertainty due to the single top 
production cross-section is shown in Table IV. Its con- 
tribution to tt and background templates is calculated 
by varying the theoretical cross-section by ±12%. 

9.7. Luminosity 

An uncertainty of 20% on the integrated luminosity 
is assumed for the early LHC data-taking period. The 
luminosity error is directly translated into the cross- 
section error and results in ±22% for the ti cross- 
section measurement [4] . 

Tables IV summarizes the contributions from each 
source of the considered systematic uncertainties on 
the selection efficiency as well as due to changes in 
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the likelihood shapes. One can sec that the largest 
contribution comes from the JES. The total system- 
atic uncertainty is calculated as a quadrature sum of 
uncertainties due to various sources. 



Table IV Summary of systematic uncertainties 



Source of systematics 


jets channel 
(%) 


e± jets channel 
(%) 


Trigger 


±1.0 


±1.0 


Lepton ID 


±1.0 


±1.0 


JES (5% up) 


±9.7 


±10 


JES (5% down) 


-11.1 


-10 


JES (10% up) 


±19.4 


±20 


JES (10% down) 


-20.8 


-21.6 


MC Generator 


±5.0 


±3.5 


ISR 


±7.7 


±10.2 


FSR 


-7.9 


-9.1 


PDF 


±1.2 


±1.6 


Single top 


±0.5 


±0.7 


Luminosity 


±22 


±22 



10. Conclusions 

The tt cross-section measurement using the likeli- 
hood method has been discussed. During early LHC 
operation the ATLAS detector will not be well under- 
stood. No b-tagging is used in this analysis. We have 
chosen likelihood variables that do not depend on jet 
energy scale and missing energy so that we can use 
them during early LHC operation and we can mini- 
mize systematic uncertainties. As a result, the domi- 
nant systematic error (jet energy scale) is reduced by 
a factor of two compared to a cut-based analysis [4]. 
The ti cross-section measurement at 50 pb~^ using the 
likelihood method can be measured with the following 
accuracy: 



fjets channel (for ± 5% JES): 

= ±5.3%(stat.) t\ll%{syst.) ± 22%{lumi) 



• ju±jets channel (for ± 10% JES): 



<^tt 



±5.3%{stat.) I24:i%(sj/si.) ± 22%{lumi 



• e±jets channel (for ± 5% JES): 



^ = ±6.0%istat.) tlli%isyst.) ± 22%{h 
<^tt 



• e±jets channel (for ± 10% JES): 



^ = ±6.0%(siai.) tlti%{syst.) ± 22%{lumi 



7) 



9) 
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